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Enhanced Stability of Layered Phases in Parallel Hard-Spherocylinders Due to 

Addition of Hard-Spheres 
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The Complex Fluid Group, Martin Fisher School of Physics, Brandeis University, Waltham MA 02254 * FOM-Institute for 
Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands 

(February 1, 2008) 

There is increasing evidence that entropy can induce microphase separation in binary fluid mixtures 
interacting through hard particle potentials. One such phase consists of alternating two dimensional 
liquid-like layers of rods and spheres. We study the transition from a uniform miscible state to this 
ordered state using computer simulations and compare results to experiments and theory. We 
conclude that (1) there is stable entropy driven microphase separation in mixtures of parallel rods 
and spheres, (2) adding spheres smaller then the rod length decreases the total volume fraction 
needed for the formation of a layered phase, therefore small spheres effectively stabilize the layered 
phase; the opposite is true for large spheres and (3) the degree of this stabilization increases with 
increasing rod length. 
PACS numbers: 64.70.Md, 64.75+g, 61.30.Cz 
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I. INTRODUCTION 

In hard particle fluids all allowed configurations have 
the same energy and therefore it is the number of states, 
or equivalently the entropy of a system that determines 
the equilibrium phase. Examples of well known phase 
transitions where the formation of ordered structures are 
driven solely by an increase in entropy are the liquid to 
crystal transition in hard spheres , the isotropic to nc- 
matic p| and the nematic to smectic transition in hard- 
rods |3||4|]. Because of their high degree of monodisper- 
sity, and because of the dominant role of steric repulsion 
in the pair-potential, colloidal suspensions of polystyrene 
latex and rod-like viruses have often been used as exper- 
imental model systems for the study of entropy induced 
ordering in hard-sphere J5|-[7| and hard-rod systems, re- 
spectively H|). 

A natural extension of the above work is to the phase 
behavior of mixtures, with a number of recent exper- 
imental and theoretical studies focusing on the phase 
behavior of binary mixtures of hard-spheres [JlO] p"9[ . 
We have recently begun work on less studied sys- 
tems that closely approximate hard-rod/hard-sphere and 
hard-rod/polymer mixtures [^o|-^6|. As a model for 
hard-rods we used either fd or TMV virus, as hard- 
spheres we used polystyrene latex, and as polymers 
we used poly(ethylene-oxide) with varying molecular 
weights |27],|2^]. The part of the phase diagram explored 
consisted of pure rods in either the isotropic, nematic, or 
smectic phase to which a small volume fraction of spheres 
or polymers was added. Remarkably, besides the ex- 
pected uniform mixtures and bulk dcmixing, we also ob- 
served a variety of microphases for a wide range of sphere 
sizes and concentrations [p7| . In microphase separation 
the system starts separating into liquid-like regions that 
are rich in either spheres or rods. However, unlike bulk 
demixing where rod and sphere rich regions grow until 
reaching macroscopic dimensions, in microphase demix- 



ing these liquid-like regions increase only to a critical 
size after which they order into well defined three dimen- 
sional equilibrium structures. One of the micro-separated 
phases observed, named the lamellar microphase, consists 
of alternating two-dimensional liquid-like layers of rods 
and spheres and is the subject of theoretical analysis in 
this paper. 

In this paper we use the second virial approximation 
first studied by Koda et. al. to examine the influence 
of molecular parameters such as shape and size, on the 
phase behavior of rod/sphere mixtures. As the second- 
virial theory is approximate in nature, we validate the 
theoretical predictions by comparing them with cither 
computer simulations or experimental results. The re- 
mainder of this paper is organized as follows: In section 
[n| we formulate the second virial approximation for the 
rod/sphere mixture. The general features of the phase 
diagram are discussed and a physical picture of the fac- 
tors responsible for the enhanced stability of the layered 
phase due to the presence of spheres is presented. In Sec- 
tion III the influence of varying the spherocylinder length 



on the phase behavior of spherocylinder/sphere mixtures 
is studied using computer simulations and the results are 
compared to theoretical predictions. Section |IV| < 



examines 



how changes of the sphere diameter influence the phase 
behavior of spherocylinder/sphere mixtures. Finally in 
Section ^ we present our conclusions. 



II. GENERAL FEATURES OF A PHASE 
DIAGRAM OF A SPHEROCYLINDER-SPHERE 
MIXTURE 

Although the equilibrium phases of all hard particle 
fluids are determined by maximizing the entropy, order- 
ing transitions are still possible because the expression 
for the total entropy, or equivalently free energy, splits 
into two parts. The ideal contribution to the entropy 



is of the form p In /?, where p is the density distribution 
function. This contribution to the entropy attains a max- 
imum for a uniform density distribution and therefore al- 
ways suppresses transitions from uniform to modulated 
phases. In contrast, excluded volume entropy sometimes 
increases with increasing order and therefore drives the 
system towards a modulated phase. In this paper we use 
a highly simplified second virial approximation to calcu- 
late the excluded volume entropy. 

The equilibrium phase in a spherocylinder/sphere mix- 
ture is determined by four parameters: length over di- 
ameter of a spherocylinder (L/D sc ), diameter of sphero- 
cylinder over diameter of sphere (D sc /D sp ), total volume 
fraction of spheres and spherocylinders (77) and partial 
volume fraction of spheres (p S p)- To help us in interpre- 
tation of our results we first define the slope 

p S p^o p sp 

where rj(p sp ) is the total volume fraction of the rod- 
sphere mixture at the layering transition after the in- 
troduction of spheres at partial volume fraction p sp . A 
positive value of r implies that adding a second compo- 
nent stabilizes the nematic phase by displacing the smec- 
tic transition to higher densities. For the case when both 
components are spherocylinders of different lengths but 
with the same diameter, slope r is positive if the ratio of 
lengths is less then approximately 7 (3^j3l[]. In the same 
manner, negative values of r imply that the second com- 
ponent stabilizes the smectic phase. There are predic- 
tions of a negative value of r in a bidisperse rod mixture 
when the ratio of rod lengths is large enough |3(| , or when 
added rods have a larger diameter [[32]]. In this section 
we focus on the phase behavior of the spherocyinder- 
sphere mixture for the specific microscopic parameters 
L/D sc — 20 and D sc /D sp — 1. We present a physical 
picture of excluded volume effects that are responsible 
for the enhanced stability of the lamellar phase. In the 
next two sections we extend our study on how changes 
in the molecular parameters L/ D sc and D sc /D sp modify 
the phase behavior and in particular, their influence on 
the magnitude and sign of the slope r. 

A. Second virial approximation 

The second virial approximation for a mixture of per- 
fectly aligned spherocylinders and spheres of equal di- 
ameter was proposed by Koda, Numajiri and Ikeda ]29| ] 
and is generalized for arbitrary L/D sc and D sc /D sp in 
the appendix. It was previously shown that the second 
virial approximation described qualitatively the forma- 
tion and various features of the smectic phase of hard 
rods P,|33|-|35| . Here we study how the addition of spheres 
perturbs the formation of the smectic phase. Since the 
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FIG. 1. Stability boundaries for a mixture of perfectly 
aligned spherocylinders (L/D 3C — 20) and spheres with equal 
diameter (D S c/D sp — 1). The full line indicates the theo- 
retical prediction of the volume fraction at which the system 
becomes unstable towards lamellar fluctuations. The dashed 
line indicates instability towards demixing into two macro- 
scopically distinct phases. Squares are results of computer 
simulations at which the layering transition is observed. The- 
oretically, the periodicity associated with a one-dimensional 
lamellar instability continuously grows and diverges as the 
system completely phase separates. Illustrations of the misci- 
ble and lamellar phases are shown in Fig. 5a and 5b, respec- 
tively. 

sphere volume fraction is very low we expect that the sec- 
ond virial approximation is still qualitatively correct for 
these mixtures. We consider a sinusoidal perturbation 
from the uniform density for both spherocylinders and 
spheres. From equations ([&]) and (|^) in the appendix we 
obtain the free energy difference between the uniformly 
mixed and layered state in a spherocylinder/sphere mix- 
ture: 

5F = a\ (su-2^Su+ (^) 2 ^ =0 (2) 

The phase diagram obtained within this approximation 
for microscopic parameters L/ D sc — 20 and D sc /D sp = 1 
is shown in Fig. [l| From the phase diagram we see that 
the first prediction of the model is that spheres, upon 
addition to a smectic phase, will preferentially occupy 
space between smectic layers and therefore create a sta- 
ble micro-separated lamellar phase. The second predic- 
tion is that the total volume fraction at which the system 
undergoes a transition from a uniform miscible state to a 
layered lamellar state is lowered by increasing the partial 
volume fraction of spheres. This implies that the slope t 
is negative for this particular spherocylinder/sphere mix- 
ture and we conclude that in this case spheres enhance 
the layering transition. 

We can assign a simple physical origin to every term 
given in Eq. (||) above and Eq.(^) of the Appendix. 
The parts of the spherocylinder-spherocylinder interac- 
tion term S22 and sphere-sphere term Sn that scale as r\ 



are due to the ideal (id) contribution to the free energy, 
also known as the entropy of mixing and are denoted as 
S22 and Slf, respectively. The terms having a rj 2 de- 
pendence in S22 , <Si2, are due to the spherocylinder- 
spherocylinder, spherocylinder-sphere and sphere-sphere 
excluded volume (ex) interaction, respectively and are 
denoted as , Sf% and Sff. Since the instability is de- 
fined as 8F(rj c ^ k c ) — 0, at a critical density r\ c and at a 
critical wavevector k c all individual contributions to the 
free energy difference in Eq. (||) must add up to zero. In 
Fig. H we show the value of all terms with distinct physi- 
cal origins at the instability density rj c and wavevector k c 
as a function of partial volume fraction of spheres. Since 
from our analysis we cannot determine the absolute am- 
plitude of a\ we only plot the ratios of all free energy 
components to the absolute value of the spherocylinder- 
spherocylinder excluded volume | |- If we set the 
partial volume fraction of spheres to zero (p sp = 0) in 
Eq.(|J) we obtain an equation whose solution indicates 
the nematic-smectic stability limit in a pure suspension 
of aligned spherocylinders (3^] . For these conditions the 
only two nonzero components of free energy are S^, 
which is negative and therefore drives the transition and 
$22, which is positive and therefore suppresses the tran- 
sition. As we start increasing the partial sphere volume 
fraction p sp , the spherocylinder-sphere free volume term 
S'ff rapidly assumes large negative values as evidenced 
by the rapidly decreasing ratio of Sff / | S22 I • This im- 
plies that layering the mixture significantly decreases the 
excluded volume that is due to the spherocylinder-sphere 
interaction. 

We can use the information gained from the second 
virial approximation to obtain a clear physical picture 
of excluded volume effects in spherocylinder/sphere mix- 
tures and explain the enhanced stability of the lamellar 
phase. Taking any single spherocylinder in a uniform 
spherocylinder/sphere mixture and replacing it by two 
spheres will leave the value of excluded volume virtu- 
ally unchanged. The reason for this lies in the fact that 
the volume excluded to the spherocylinder due to the 
presence of a sphere with equal diameter, under the con- 
straint of uniform packing, is a spherocylinder with di- 
ameter 2D SC and length (L + 2D SC ) where L and D sc are 
defined in Fig. ^| However, the excluded volume between 
any two spherocylinders with large L/D sc is only about 
twice this value as illustrated in Fig. [jj. Although replac- 
ing spherocylinders by spheres in such a manner leaves 
the excluded volume almost unchanged, it does signifi- 
cantly decrease the total volume fraction of the mixture 
since the volume of two spheres is much smaller then the 
volume of a spherocylinder with large L/D sc . Therefore 
in the spherocylinder/sphere mixture we encounter ex- 
cluded volume problems similar to those found in a pure 
spherocylinder solution, but at a lower total volume frac- 
tion. As in pure spherocylinders, the system reduces the 
excluded volume by undergoing a transition to a layered 
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FIG. 2. Term-by-term dependence of the free energy dif- 
ference between the miscible and lamellar phases (Eq. ^) 
as a function of the partial volume fraction of spheres for 
L/D ac = 20 and D sc /D 3p = 1. The Slf and S22 terms are 
the sphere and spherocylinder ideal contributions to the total 
free energy difference between the layered and uniform states. 
Sif, S12 and 5*22 are excluded volume contributions to the 
free energy due to sphere-sphere, spherocylinder-sphere and 
spherocylinder-spherocylinder interactions respectively. Since 
from our analysis we cannot determine the amplitude in Eq. 
^| we plot amplitude independent ratios of each of five com- 
ponents of the free energy with different origins to the sphe- 
rocylinder-spherocylinder excluded volume interactions. The 
stability condition is that SF = 0, so for any value of partial 
volume fraction of spheres p ap the sum of the five contribu- 
tions to SF is zero. SF of the ideal terms are positive, hence 
they stabilize the uniform, miscible nematic state, while the 
free volume terms are negative, favoring the lamellar state. 
The excluded volume sphere-sphere term (5*if)is negligible 
and the spherocylinder-sphere (S12 ) term dominates the tran- 
sition. 
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FIG. 3. a) Volume excluded to the center of mass of a 
spherocylinder (sc) due to the presence of a sphere (sp) is 
indicated by light shading b) Volume excluded to the center 
of the mass of a second spherocylinder due to the presence 
of the first. Replacing a spherocylinder by a sphere decreases 
the excluded volume by approximately a factor of two, but 
it decreases the total volume fraction much more since the 
volume of a spherocylinder with large L/D ac is greater than 
the volume of a sphere with diameter D sc . The comparatively 
large excluded volume between a sphere and a spherocylinder 
is the reason for the enhanced formation of the lamellar phase 



phase. The excluded volume is reduced in the lamel- 
lar state because a periodic density distribution forces 
spheres and spherocylinder into alternate layers thus de- 
creasing the probability of the very unfavorable sphere- 
spherocylinder contacts as illustrated in Fig. |^. This 
explains the large decrease in the value of the Sjjf term 
at the lamellar transition that we observe in the second 
virial theory This term is responsible for the enhanced 
stability of the lamellar phase in a sphere/spherocylindcr 
mixture. In conclusion, it is the inability to efficiently 
pack a uniform mixture of sphcrocylinders and spheres, 
as reflected in the large spherocylinder/sphere excluded 
volume term, that destabilizes the nematic phase and en- 
hances the formation of a layered phase. 

An alternate way to think about the formation of a 
layered phase is to focus on the effects of spherocylinder 
ends . The nematic phase in our simplified model 
is characterized by random distribution of spherocylin- 
ders along their axial and radial direction as illustrated 
in Fig. H This end effect is responsible for the for- 
mation of the smectic phase, which is characterized by 
a periodic density distribution. In similar fashion, in- 
troducing a sphere into the nematic phase will have the 
same effect on the surrounding spherocylinders as an- 
other spherocylinder end. Therefore adding spheres very 
effectively increases the density of "spherocylinder ends" 
and decreases the total volume fraction. To resolve the 
difficulties in efficient packing due to these extra "sphero- 
cylinder ends" , the mixture layers at a lower total volume 
fraction. 



B. Monte Carlo Simulation 

In the previous section we discussed two predictions 
of the second virial theory for a spherocylinder/sphere 
mixture with L/D sc = 20 and D sp /D sc = 1; the exis- 
tence of the lamellar phase and the enhanced stability 
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FIG. 4. a) A schematic illustration of excluded volume 
effects in a nematic phase in a spherocylinder/sphere mix- 
ture. In the nematic or miscible phase each sphere creates 
a large excluded volume around it, indicated by gray areas, 
that is inaccessible to spherocylinders. b) When the system 
undergoes a transition to a layered phase, the large excluded 
sphere-spherocylinder volume vanishes since the probability 
distribution severely limits the number of ways that spheres 
are allowed to approach spherocylinders. 



of the lamellar phase when compared to a smectic phase 
of pure spherocylinders. Our results are in agreement 
with previous studies of Koda et. al. (^9|. However, the 
second virial approximation is highly approximate and 
there is reasonable concern about the influence of higher 
terms on the topology of the phase diagram. To sup- 
port their conclusions Koda et. al. performed computer 
simulations, which indicated the existence of an lamel- 
lar phase |37||29| ]. Still, the question of whether spheres 
simply fill the voids between layers in an already formed 
smectic phase, or actually induce layering at lower total 
volume fraction was not addressed. In this section, us- 
ing Monte Carlo simulations we address the question of 
the influence of adding spheres on the phase behavior of 
spherocylinders by determining the slope r in Eq. [I] in a 
mixture of spherocylinders and spheres with parameters 
L/D sc = 20 and D sp /D sc = 1 

A Monte Carlo simulation of a mixture of hard- 
spheres and perfectly aligned hard-spherocylinders was 
performed at constant pressure and number of parti- 
cles Q . Most simulations contained 392 spherocylinders 
and a variable number of spheres. To check for finite size 
effects we also ran simulations with 784 spherocylinders, 
but saw no significant difference in the results obtained. 
In one sweep, pressure was increased from a dilute ho- 
mogeneous mixture up to a well ordered, dense smectic 
or lamellar phase. At each value of the pressure, the 
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FIG. 5. a) A schematic example of a typical configura- 
tion of spherocylinders in a dense nematic phase. Since the 
nematic phase is characterized by a uniform density distri- 
bution this results in inefficient packing and large excluded 
volume between spherocylinders both along their radial and 
axial directions. This large and unfavorable excluded volume 
is indicated by lightly shadowed areas, b) An illustration of a 
typical configuration of spherocylinders in a columnar phase 
where the excluded volume between spherocylinders is lower 
compared to the nematic phase at the same density, and the 
ideal part of free energy is higher. In a columnar phase the 
spherocylinders are forced into registry as one spherocylinder 
occupies space right above or below another one. Therefore 
the columnar phase is characterized by two dimensional order 
in the plane perpendicular to the spherocylinder's long axis 
and one dimensional disorder parallel to the long axis, c) A 
representative configuration of spherocylinders in a smectic 
phase, which is characterized by one-dimensional order along 
the spherocylinder's long axis and two-dimensional disorder 
in the perpendicular directions. Both theory and experiment 
indicate that the columnar p hase i s always metastable with 
respect to the smectic phase 



E 2? 

03 CD 




smectic 



nematic 



0.25 0.30 0.35 0.40 0.45 0.50 0.55 
total volume fraction 
FIG. 6. Smectic order parameter obtained from Monte 
Carlo simulations is plotted against the total volume fraction 
for spherocylinders with L/D ac — 20. From right to left, the 
partial volume fraction of spheres (p sp ) increases from 0% 
to 2.4% as indicated by the legend. The phase diagram is 
reconstructed from this data by defining a phase as layered 
when the spherocylinder order parameter reaches a value of 
0.3. 



density of spheres and spherocylinders and their corre- 
sponding smectic order parameter were measured after 
the system was allowed to equilibrate. Identical results 
were obtained when the pressure was slowly decreased 
from a initially dense phase composed of alternating lay- 
ers of spherocylinders and spheres to a dilute homoge- 
neous mixture. 

Besides lamellar transitions there is a possible demix- 
ing transition where spherocylinders and spheres phase 
separate into macroscopically distinct phases. However, 
once a layered phase is formed the exchange of spheres be- 
tween layers drops to a negligible amount, leaving open 
the possibility that system would undergo a demixing 
transition, but is stuck in a lamellar phase, which is only 
a metastable state. To find out the location of the demix- 
ing transition it is necessary to measure the chemical po- 
tential of both spherocylinders and spheres in a sphero- 
cylinder/sphere mixture (39). This possibility was not 
examined in this work, primarily because we are only 
interested in how low concentrations of spheres perturb 
the formation of the layered phase. Therefore it is rea- 
sonable to expect that at a very low volume fraction of 
spheres, the lamellar transition is going to be more stable 
than the demixing transitions as predicted by the second 
virial theory. 

A plot of the smectic order parameter for spherocylin- 
ders with L/D sc — 20 as a function of increasing total 
density for different partial volume fractions of spheres is 
shown in Fig. ^. As the system approaches a certain crit- 
ical density we observe a rapid non-linear increase in the 
smectic order parameter that we interpret as a signature 
of the nematic to smectic phase transition. This critical 
density shifts to lower values of the total volume fraction 
as the partial volume fraction of spheres is increased. To 
reconstruct a phase diagram from the above data we de- 
fine a phase as layered when its smectic order parameter 



reaches a value of 0.3 [Icj ]. For a pure spherocylinder 
suspension this value yields good agreement with previ- 
ous studies of the volume fraction of the nematic-smectic 
phase transition |nj. Since we are mostly interested 
in the qualitative behavior of a spherocylinder/sphere 
mixture this method should suffice our purposes. Us- 
ing this phcnomcnological rule, the phase diagram for 
a mixture of spherocylinders and spheres (L/D sc = 20, 
D sc /D sp — 1) is reconstructed and compared to the sec- 
ond virial theory in Fig. [I]. An immediate conclusion 
drawn from Fig. [I] is that adding spheres to aligned sphe- 
rocylinders enhances the stability of the lamellar phase, 
which is indicated by the negative value of slope r, in 
agreement with the prediction of the second virial ap- 
proximation. 



III. THE EFFECTS OF SPHEROCYLINDER 
LENGTH ON THE PHASE DIAGRAM 

Next we proceed to investigate the influence of vary- 
ing the spherocylinder length on the magnitude of slope 
t. The predictions of the second virial theory for the 
nematic-lamellar instability are shown in Fig. ^|a. The 
second virial theory clearly predicts increasing stability of 
the lamellar phase with increasing length of spherocylin- 
der. To verify this prediction we repeated Monte Carlo 
simulations for spherocylinders with different L / D sc and 
used the same rule as before to identify the volume frac- 
tion of the nematic-lamellar transition. The simulation 
results for the location of the nematic to layered tran- 
sition are shown in Fig. [?]b. We can conclude that 
our simulations confirm predictions of the second virial 
model and that the length of the spherocylinder is an 
important parameter in forming the lamellar phase, with 
longer spherocylinders showing an increasing tendency to 
form a layered phase at a lower volume fraction of added 
spheres. 

Using the physical picture of the excluded volume ef- 
fects developed in the previous section provides a natural 
explanation for our simulation results in Fig. |^. With in- 
creasing spherocylinder length the excluded volume due 
to the spherocylinder-sphere interaction grows propor- 
tionally to the spherocylinder length and consequently 
the value of the Sjf term increases in magnitude. As we 
have seen before, the larger the Sff term, the more likely 
it is for the system to form a layered phase. 

It is interesting to consider the limit of spherocylinders 
with infinite aspect ratio. In the density regime of the 
nematic-smectic transition, this model can be mapped 
onto a system with skewed cylinders with an aspect ra- 
tio close to one. The nematic-smectic transition in this 
model has been studied numerically |^,^|. If we con- 
sider the addition of spheres to this system, then the same 
affine transformation that maps the infinite spherocylin- 
ders onto squat, skewed spherocylinders, will map the 
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FIG. 7. a) Prediction from the second virial theory for the 
total volume fraction 77 of the lamellar instability as a func- 
tion of sphere partial volume fraction (p sp ) for spherocylinders 
with different L/D 3C ratios. The diameter of spherocylinders 
is kept constant and is equal to the diameter of the spheres, 
b) Results from Monte Carlo simulations for the lamellar in- 
stability of spherocylinders as a function of partial volume 
fraction of spheres for same conditions as in Fig. ^a. The 
volume fraction at the phase transition was defined as having 
a smectic order parameter of spherocylinders equal to 0.3 



spheres onto infinitely thin, parallel disks. As the disks 
are infinitely thin, they do not interact with each other 
but only with the cylinders. Inside the nematic phase, 
most volume is excluded for these disks. However, in the 
smectic phase, there is ample space for the disks between 
the layers. In fact, the stronger the layering, the larger 
the accessible volume. Hence in this limit, the addition 
of spheres will strongly stabilize the smectic phase. 

IV. THE EFFECTS OF SPHERE DIAMETER ON 
THE PHASE DIAGRAM 

In this section we investigate the influence of sphere 
diameter on the value of slope r. Fig. || shows the 
prediction of the second virial theory for the depen- 
dence of slope t on the ratio of spherocylinder to sphere 
diameter (D sc / D sp ) for spherocylinders with different 
L/D sc . In section A we examine the phase behavior of 
sphere/spherocylinder mixtures when the sphere diame- 
ter is smaller then spherocylinder diameter and in section 
B we examine the other case when the sphere diameter 
is larger then the spherocylinder diameter. In our model 
the presence of the spheres cannot alter the orientational 
distribution function of spherocylinders, which are al- 
ways perfectly parallel to each other. It is reasonable 
to expect that this assumption holds for spheres smaller 
then the spherocylinder length, but as a sphere becomes 
larger then the spherocylinder length, long wavelength 
elastic effects start to dominate the behavior of the sys- 
tem and hard spherocylinders will tend to align parallel 
to the surface of the sphere Q]. Therefore in Fig. || we 
plot the values of slope r only for those values of D sc / D sp 
for which our assumptions are at least qualitatively cor- 
rect. As we increase the sphere size beyond this limit our 
model describes a highly artificial system of large spheres 
and parallel spherocylinders. In this regime we observe 
oscillations in the value of slope r similar to what is ob- 
served in binary mixtures of parallel spherocylinders |30| . 

A. Sphere diameter smaller than spherocylinder 
diameter 

In the regime where D sc / D sp > 1 (for spherocylinders 
of any L/D sc ), decreasing the sphere size increases the 
stability of the lamellar phase as indicated by the in- 
creasing negative value of slope r seen in the right hand 
side of Fig. ||. This prediction of the theory has a sim- 
ple explanation in our picture of excluded volume in a 
sphere/spherocylinder mixture. If we halve the sphere 
radius D sp , while keeping constant the volume fraction 
of spheres, we increase the number of spheres eight times. 
At the same time, the result of reducing the sphere size 
is to decrease the excluded volume of the spherocylinder- 
sphere interaction. However, the eightfold increase in the 
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FIG. 8. Theoretical prediction for the stability criterium 
of the lamellar phase r in Eq. ^ as a function of spherocylin- 
der (sc) to sphere (sp) diameter ratio for four spherocylinders 
with different L/D sc . The negative value of slope r indicates 
that spheres of that particular size enhance the layering tran- 
sition. Larger negative values of r implies the formation of 
the lamellar phase at a lower total volume fraction. 

number of spherocylinder-sphere interactions more then 
compensates for the decrease in excluded volume between 
the sphere and spherocylinder and consequently the mag- 
nitude of Sff increases with decreasing sphere diameter. 
This leads to the increased stability of the layered phase 
with decreasing sphere size. 

It becomes difficult to verify this prediction using com- 
puter simulations. As the sphere size decreases at con- 
stant total volume fraction r], the number of particles in 
a simulation rapidly reaches the order of thousands re- 
quiring simulation times that are prohibitively long. As 
the ratio of spherocylinder to sphere diameter (D sc /D sp ) 
was varied within the accessible range between 0.5 to 2 
we did not observe any changes in the value of slope r 
that were larger than our measurement error. Larger 
and longer simulations are needed for a careful analysis 
of spherocylinder/sphere mixtures with extreme values 
of the ratio D sc /D sp . 

B. Sphere diameter larger than spherocylinder 
diameter 

For spherocylinders with small L/D sc , Fig. || shows 
that the magnitude of slope r uniformly decreases with 
increasing sphere size. Eventually the slope r changes 
sign and becomes positive, implying that large spheres 
stabilize the nematic and not the smectic phase. The 
phase diagram under conditions where slope r is positive 
is shown in Fig. |^. The wavevector associated with the 
layering transition, indicated with a solid line in Fig. |^, 
remains at an almost constant value. Another impor- 
tant point is that the amplitude ratio in Eq. (j^) is pos- 
itive. This means that the periodic density modulations 
of the spherocylinders and spheres are in phase, which 
implies that spheres no longer go into the gap between 
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FIG. 9. Stability diagram of a mixture of spherocylinders 
(L/D ac = 10) and large spheres D 3C /D sp = 0.15. Unlike a 
mixture of small spheres and spherocylinders (Fig. [j]), intro- 
ducing large spheres displaces the layering transition to higher 
total volume fractions indicating a positive value of slope r. 
The structure of the layered phase is also different, with large 
spheres fitting in the smectic layer rather than into the smec- 
tic gap. The smectic periodicity associated with the layering 
transition does not change significantly until the concentra- 
tion of spheres is high enough for the system to demix. Then 
the smectic wavevector discontinuously jumps to zero. 

two spherocylinder layers, but rather fit into the sphcro- 
cy Under layer. However, as the partial volume fraction 
of spheres (p sp ) is increased further we observe a dis- 
continuous jump in the wavevector to zero value. This 
implies that there is a discontinuous change from a layer- 
ing to a demixing transition. As the demixing transition 
is reached there is also a change in sign of the amplitude 
ratio, which becomes negative and the spherocylinders 
and spheres bulk separate. In contrast, the phase dia- 
gram for mixtures of small spheres and spherocylinders 
shown in Fig. [l] looks quite different. The amplitude ra- 
tio for this case is always negative implying formation of 
the lamellar phase. Another contrast is that in a mix- 
ture of small spheres and spherocylinders the wavevector 
associated with the layering transition decreases in a con- 
tinuous fashion until it reaches zero value. 

We now examine the behavior of individual terms in 
Eq. (||) for a mixture of large spheres and short sphe- 
rocylinders shown in Fig. ||. Most notably, we find that 
at low volume fractions of spheres where the system un- 
dergoes the layering transition, the ratio Sff/S'If << 1- 
This implies that upon layering there is almost no reduc- 
tion of the unfavorable sphere/sperocylinder interaction 
and that the spherocylinder/spherocylinder interaction 
alone drives the formation of the layered phase. In con- 
trast, for small spheres this ratio was large and was re- 
sponsible for enhanced stability of the lamellar phase as 
was shown in Fig. |^. At a higher volume fraction of large 
spheres where the mixture directly bulk phase separates 
we find that the ratio S'ff/S'ff >> 1. This implies, as 
expected, that demixing very effectively reduces the un- 
favorable sphere/spherocylinder interactions. These re- 



sults suggest a physical picture of the excluded volume 
effect. Unlike small spheres, large spheres can not fit into 
the gap between smectic layers and consequently there is 
no way to gain free volume by undergoing the layering 
transition. As an alternative, to gain free volume the sys- 
tem bulk phase separates at the lowest volume fraction 
of spheres possible. 

While for short spherocylinders the magnitude of slope 
t uniformly decreases with increasing sphere size, longer 
spherocylinders exhibit a qualitatively different behavior. 
For a mixture of spherocylinders with L/D sc = 100 and 
spheres with D sc /D sp = 0.1 there is a pronounced in- 
crease in the stability of the lamellar phase as shown in 
Fig. By increasing the length of spherocylinders to 
even larger values, the region of increased stability of the 
lamellar phase shifts to higher values of the sphere ra- 
dius. Two conditions emerge, which when satisfied lead 
to enhanced stability of the lamellar phase. First, it is 
necessary for a sphere to fit between two smectic lay- 
ers without disturbing them. This condition is satis- 
fied when D sp /L w 0.1. The second condition is that 
D sp /D sc >> 1. It was argued before that under these 
condition large spheres are able to induce smectic corre- 
lations amongst neighboring spherocylinders p7[ , which 
in turn can enhance the formation of the lamellar phase. 

Because of the large size asymmetry it was not feasi- 
ble to carry out simulations for mixture of spherocylin- 
ders and spheres with L/D sc « 100 and D sc /D sp « 0.1. 
However, these conditions are closely approximated by 
recent experiments on rod-like fd (L = 1/xm, L/D sc w 
100) and polystyrene spheres |2^| . Therefore, we com- 
pare theoretical results of slope t for spherocylinders 
with L/D sc = 100 shown in Fig. || to these experi- 
mental results (27J. When large spheres D sp w l//m, 
(D sc /D sp w 0.01) are mixed with fd at any concentra- 
tion for which the nematic phase is stable, we observe 
no formation of the layered phase. Instead, large spheres 
phase separate into dense aggregates elongated along the 
nematic director indicating that the value of slope r is 
larger then zero. When the size of the sphere was de- 
creased to D sp = 0.1/Mn, (D sc /D sp ps 10) we observed 
a transition to a layered state at a fd concentration of 
20 mg/ml. The formation of a smectic phase in a pure 
fd suspension at the same ionic strength occurs at 65 
mg/ml. The fact that adding spheres diminishes the rod 
density by a factor of three indicates a large negative 
value of slope r. As the sphere size was further decreased 
D sp = 0.022^m, (D sc /D sp — 0.46) there was again indi- 
cation of a lamellar phase, but this time at a much higher 
concentration of rods of about 50 mg/ml. Thus, although 
small spheres still stabilize the layering transition, imply- 
ing a negative value of slope r, the magnitude of slope r 
is much less for D sc /D sp « 0.46 then for D sc /D sp m 0.1. 
These qualitative trends of the non-monotonic behavior 
of slope t with sphere size observed in experiments of 
/fi-polystyrene mixtures are very similar to the theoreti- 



cal prediction shown in Fig. 
L/D sc = 100. 



for spherocylinders with 



V. CONCLUSIONS 

In this paper we have presented the predictions of 
the second virial theory for a mixture of parallel hard- 
spherocylinders and hard-spheres undergoing one dimen- 
sional microphase separation. We have been able to ver- 
ify a number of these predictions using Monte Carlo simu- 
lations. We found that spheres induce layering, which im- 
plies a negative value of the slope r, which is the change 
in total volume fraction of the mixture at the point of 
nematic-smectic instability with respect to the partial 
volume fraction of added spheres (Eq. ||) . At the same 
time the magnitude of the slope r increases with increas- 
ing spherocylinder length. In other words, spheres at the 
same partial volume fraction stabilize layering of longer 
spherocylinders more then shorter spherocylinders. Be- 
sides this, the theory predicts an unusual non-monotonic 
behavior in slope r as a function of sphere to sphero- 
cylinder diameter. Although the physical origin of this 
effect is not clear, it is intriguing that similar qualitative 
trends are observed in experiments of mixtures of the 
spherocylinder- like fd and polystyrene spheres. However, 
in real experiments spherocylinders are free to rotate, are 
flexible, and have charge associated with them. Before 
quantitative comparisons with experiments are possible 
it will be necessary to perform simulations and formu- 
late theories that take into account these effects mostly 
ignored in this highly idealized treatment. 
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VII. APPENDIX 

A general expression for the free energy of bidisperse 
mixture at the second virial level is 

pF{p 1 ,p 2 )=Y t f d(r) Pi (r)ln( Pi (r)) - 

„ Y\ Y) / dri / dr 2 p i (r 1 )p i (r 2 )/ iJ (r 1 ,r 2 ) (3) 

i=i, 2j =i,2^ J y 



where the function /jj is the overlap function between 
two spheres, sphere and spherocylinder or two sphero- 
cylinders p9| . It attains the value of -1 if two particles 
overlap, otherwise it is equal to 0. The terms involving 
p In p represent the entropy of mixing while the terms in- 
volving fij represent the free volume entropy. Since we 
are interested in one dimensional layering we look at the 
response of the system to following density pertubation 



Spi(z) — a\ cos(k z z) 
6p2(z) = ai cos(k z z) 



(4) 



The free energy difference between the uniform and per- 
turbed state is 

5F = F(l + 5 Pl {z), 1 + S P2 (z)) - F(l, 1) = aSa (5) 

where a = (ai,a 2 ) and S is a two dimensional stability 
matrix. To find the limit of stability we have to solve the 
equation det(S) = 0. For latter convenience we define 
the following function 

L 3sin(fccr(2 + 2^)) 



S(-^-,(T,k) = 



4fc 3 

2fcercos(fccr(2 + 2-^-)) - sin(fc2a^ 
IP 



(0) 



The above expression depends only on geometrical fac- 
tors and is related to the Fourier transform of the sphe- 
rocylinder which is specified by the excluded volume be- 
tween a sphere of diameter D sp and a spherocylinder of 
length L and diameter D sc . Wavevector k is dimension- 
less because it is rescaled with the spherocylinder diam- 
eter (D sc ). The parameter a is defined as ratio of sphere 
diameter to spherocylinder diameter (a = D sp /D sc ). In 
the limit of L/D sc — > the above expression reduces to 
a Fourier transform of a sphere with unit diameter. The 
stability matrix S for a mixture of spherocylinders and 
spheres has the following form 



/ ??(! + 4(1 - Psp )r^(0,l,fc)) 
2 Ps P (l - p sp )rj 2 S{ — , l + a,k) 



\ 



L 



2p sp (l - Ps P )v S ( jj-> 1 + CT > k ) 



VPspK 



3 L 
2~5Z 



1) 



r]p 8 pS(2—,2a,k)) 



(7) 



where p sp denotes partial volume fraction of spheres and 
varies between and 1 while 77 denotes total volume frac- 
tion. Note that the terms in matrix elements S\\ and S22 
proportional to rj are due to configurational entropy while 
terms proportional to rj 2 are due to free volume entropy. 
As k — > the condition det(S) = reduces to the usual 
thermodynamic condition for the stability of the system 
against bulk phase separation. 

To reconstruct the stability diagram from the deter- 
minant we slowly increase the total volume fraction 77. 
At a certain value of total volume fraction (r) c ) the de- 
terminant of S will equal zero for a specific wavevector 
(k c ). If the wavevector k c obtained has a finite value it 
implies that system is undergoing a layering transition. 
On the other hand, the condition det(S) = when k c = 
implies complete demixing. Once we obtain values of r\ c 
and k c we can find out the ratio of amplitudes from the 
following formula 



ai _ Si 2 (r) c , k c ) 
ci2 Sn(r] c ,k c ) 



(8) 



A positive value of the amplitude ratio implies that the 
spheres and spherocylinders are in the same layer (the 
periodic modulations are in phase), while a negative value 
implies that the spheres and spherocylinders intercalate 
(the periodic modulations are out of phase). 
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